PROBING DENSITY FLUCTUATIONS USING THE FIRST RADIO 
SURVEY 



C.M.CRESS + FIRST COLLABORATION 
Columbia University 

Astronomy Dept., 538 W120th St, New York, NY10027, USA 



1. Introduction 

For decades, the clustering of luminous objects has been used as an indirect probe 
of the density fluctuations in the universe. Although the relationship between 
luminous objects and mass in the universe is not fully understood, it has been 
hoped that the observed evolution of clustering of different luminous populations 
might constrain many cosmological unknowns - the geometry of the universe, the 
nature of any non-luminous matter, the nature of density fluctuations in the early 
universe and hence the physics of the early universe. Indeed, models for structure 
formation which are based on variants of the 'cold dark matter' (CDM) model 
have received much attention largely because they have been fairly successful at 
predicting the clustering of luminous objects. 

Much information on the clustering of nearby populations has been obtained 
from optical and infra-red surveys of galaxies (and of clusters of galaxies). In 
addition, Peacock and Nicholson (PN, 1991) investigated the clustering of bright 
(S2.7 > Uy) nearby (0.01 < z < 0.1) radio galaxies. The new radio surveys provide 
an opportunity for investigating the clustering of much fainter radio sources. Prob- 
ing density fluctuations with a different tracer provides interesting information on 
the relative 'bias' of different populations and, of course, one might hope that 
the great depth and large area covered by the new surveys can be used to probe 
clustering at higher redshift and on larger scales than was previously possible. 

2. Measuring the correlation function in the FIRST Survey 

The FIRST survey has now been extended to include a total of 3000 deg 2 of sky 
yielding a catalog of about 250,000 sources with 22° < S < 42°. and 7 ,l 30 m < a < 
17 h 30 m . Details of the mapping and catalog generation can be found in Becker et 
al. (1996) and in White et al. (1997). Details of how the angular correlation function 
(CF) is measured can be found in Cress et al. (1996). Here, we use the CF estimator 
given by Landy and Szalay (LS, 1993): w{9) = [DD{6) - 2DR{6) + RR{6)}/ RR{9). 
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Figure 1. The correlation function for (i) all sources and (ii) for those sources which are identified 
with E<19 galaxies in the APM survey. 



While single objects are sometimes resolved into two or more sources in the cata- 
log, it appears that these extra sources at small separations do not contribute to 
the correlation function at larger separations (see Cress et al. 1996). Error bars 
are obtained using the 'partition bootstrap method' in which the CF is calculated 
for 10 subdivisions of the survey region and the standard deviation of these mea- 
surements at each angle is used as a measure of the error. All areas in the survey 
where the rms noise is less than 0.2 mJy were used in the analysis. CF parameters 
were determined from a straight line fit to the log-log plots. 

We also investigated the CF of FIRST sources in the catalog that have a galaxy 
with E < 19 within 3" in the APM catalog of the POSSI survey plates. 

Writing the angular CF as w(9) = A9 l_1 with measured in degrees, we find 
A = 0.002, 7 = 2.1 for the whole sample and A = 0.04, 7 = 1.8 for the APM 
matches. 

3. Inferring spatial information 

The spatial CF, £, is related to the angular CF via Limber's equation. 

fee ,-oc / dN ,\2 

•" (e) = 2 i I {■*■*) (lr - s,dxd "- (1) 

where the physical separation is given by r 2 = a 2 (t)[u 2 / F 2 (x) + x 2 9 2 }, u is the 
difference in radial distance to two sources, x is the average distance to the two 
sources, 9 is the angle between the sources in radians, a is the scale factor and 
F = 1 in a flat universe (Baugh and Efstathiou, 1993). 
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If £ has a power-law form and density perturbations evolve linearly then we can 
write £(r, z) = (r/ro) 7 (l + z)~ 2 (in comoving coordinates). We can then estimate 

using luminosity functions give in Dunlop and Peacock (DP, 1991) or Condon 
(1984). Using the LS estimator of the CF and DP's models 1 to 5 we obtain an 
r around 8/i _1 Mpc. DP's model 7 and Condon's model give slightly lower values: 
between 8/i _1 Mpc and 6/i _1 Mpc. (Note that the ro derived in Cress et al. uses the 
CF given by a different estimator - we later found the LS estimator to be more 
robust to changes in survey geometry). The ro results depend on the assumed 
evolution but the derived clustering appears fairly similar to the clustering of 
galaxies detected in the optical. 

Of course, the spatial CF is not necessarily well-represented by a pure power law 
and it is interesting to take spatial CF's predicted by CDM models and investigate 
what they imply for the clustering of objects with the redshift distribution of 
FIRST sources. CDM makes predictions for how the mass clustering evolves but 
one must also consider the possiblity that the luminous sources are 'biased' (Kaiser 
1987) relative to the mass distribution. 

In CDM models one starts with a scale- free primordial power spectrum, P(k) oc 
k n (where k is the comoving wavenumber). This P(k) is then 'processed' according 
a transfer function for CDM. At the epoch at which structure begins to grow 
P(k) = P k n T 2 (k), where T(k) for CDM is given by Bond & Efstathiou (1984) 
and Pq is the normalisation. We normalise to COBE 4-year data using Liddle et 
al.'s (1996) formula which is valid for spatially flat models. 

In the linear regime, the time evolution of the power spectrum can be calculated 
analytically and depends only on the expansion of the Universe. It is given by 
P(k, z) = P{k, z = 0) x G 2 (fl (z), n A (z))/(l + z) 2 where the growth factor, G, is 
given in Carroll, Press & Turner (1992). Note that G(fl = 1, Q A = 0) = 1. 

In the highly non-linear regime, the CF obtained from a scale-free primordial 
spectrum obeys a simple scaling relation (Groth & Peebles 1977). One can inter- 
polate between the linear and highly non-linear regimes to produce semi-analytic 
models for clustering in the quasi-linear regime which can then be accurately fit to 
results from A^-body simulations (Hamilton 1991; Jain, Mo, & White 1995; Pea- 
cock & Dodds 1996). Following the discussion of PD, we define the dimensionlcss 
power spectrum as A 2 (k,z) = (27r 2 ) _1 fc 3 P(fc, z). The non-linear dimensionlcss 
spectrum is then given by A 2 NL {k^L, z) — fpfL[A 2 (k 1 z)], where the linear and 
nonlinear scales are related by k = [1 + L (k^ l)]^ 1 ^ 3 kjq l ■ The form of Jnl for 
various values of f^o is given in PD. 

Once we have an expression for the power spectrum (including contributions 
from non-linear evolution) at all epochs, we can Fourier transform this and sub- 
stitute into Limber's equation to obtain an estimate of the predicted angular CF. 

Figure 2(a) shows how the predictions for w(9) vary for different estimates 
of the redshift distribution and how they change when non-linear evolution is 
considered. There appears to be a fair amount of uncertainty in the predicted 
result, depending on which model is chosen. It should be noted, however, 
that Condon's luminosity function is more carefully constructed to fit the low- z 
population and it is these sources that one expects to be contributing most to the 
clustering signal (see below). In addition, Peacock has indicated that DP's model 
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Figure 2. CDM predictions for w(6). (a) shows the effect of varying the assumed redshift 
distribution and the effect of including non-linear evolution in an fio = 1, Ho = 50 universe (bold 
curves include non-linear evolution), (b) shows the effect of varying the cosmological parameters 
using the DP7 redshift model (CF amplitude increases with increasing Ho). 

7 is the best model to use. If one limits oneself to only these two models then 
uncertainties in the redshift distribution do not contribute very significantly to 
uncertainties in the predicitions. The effect of non- linear evolution is significant 
on scales less than ~ 30'. 

4. Discussion 

It is useful to investigate how one expects the contribution to the clustering sig- 
nal to vary with redshift. Since Condon's model and DP7 include a large local 
population, these two models predict a large contribution to the clustering signal 
from 'local' sources: ~ 60% of the signal at 30' is expected to come from sources 
with z < 0.1 and this fraction increases as one goes to larger angles. A scenario 
where the very nearby population is made up of a mix of spirals and ellipticals 
(starbursting galaxies and AGN) and the slightly more distant population begins 
to be dominated by AGN is consistent with a clustering amplitude fairly close to 
that of 'normal' galaxies plus a slope more similar to that observed for elliptical 
galaxies (7 ~ 2.1 for ellipticals). 

Figure 2(b) shows the effect of varying cosmological parameters on the pre- 
dicted CF's but, given the above comments, it is somewhat misleading. It appears 
that an fio = 1 model with a high Hubble constant fits the data best, but this 
assumes that all the sources have similar clustering strengths. On large angu- 
lar scales (larger than ~ 30'), however, the CF is probing nearby structures and 
seems to be indicating a bias, b, for these sources of a little more than 1 (where 
iradio = b 2 ^ mass ) . On the smaller scales where one expects to be more sensitive to 
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higher- z clustering, a larger bias is required if models with popular values of Ho 
are to fit. This is consistent with the difference in clustering strengths measured by 
PN for local AGN (10/i _1 Mpc) and that measured for optically observed galaxies. 

The angular CF of radio sources associated with E > 19 galaxy matches is 
significantly higher than what one would expect for a normal sample of E > 19 
APM galaxies. This can be explained if the APM galaxies with radio counterparts 
are, on average, closer than all APM galaxies with E > 19. Ultimately, one would 
like to investigate clustering at high- z, but this is difficult given that the angular 
CF measured here appears rather sensitive to the clustering of a nearby popula- 
tion of sources which probably does not have the same clustering strength as the 
population at higher-z. In a preliminary attempt to isolate the clustering of the 
AGN-dominatcd sample at higher-z, we determined the CF for all sources without 
galaxy counterparts. The r Q 's inferred from this measurement (with a linear evolu- 
tion assumption) ranged from 10ft- _1 Mpc to 15/i _1 Mpc depending on the redshift 
model. If the redshift distribution were better determined in the moderate-z range 
(beyond what one expects for APM galaxies, but not so far that the projected 
clustering becomes negligible) , we could compare the clustering in this range with 
PN's measurement for local AGN to investigate the evolution of clustering in a 
more homogenous population. 

5. Weak Lensing 

Since it is not clear how light traces mass, it is desirable to find methods of probing 
mass fluctuations directly. In optical images, the weak gravitational lensing (co- 
herent distortion) of background sources by foreground matter distributions has 
been used to probe cluster-scale masses and the use of similar techniques to probe 
larger mass scales has been discussed (see Kaiser 1996 and references therein). 
Here, we outline a preliminary search for a weak lensing signal in FIRST. A de- 
tailed discussion will be given in Refregier et al. (1997). 

We start by defining the ellipticities of the sources. If a and b are the major 
and minor axes of a source and a is its position angle measured relative to some 
arbitrary set of axes defined on the sky, then stretching along these axes can be 
measured by e+ = e cos(2a) and stretching at 45° to these axes is measured by e x = 
ecos(2a) [where e = (a 2 — b 2 ) /(a 2 + b 2 )]. It is helpful to define ellipticity correlation 
functions which are independent of the coordinate system. To do so, we define e r + 
and e r x measured with respect to axes which are parallel and perpendicular to the 
line connecting 2 points which are being correlated. If <j) is the angle between the sky 
axes and the axes defined by the two points then one gets e T + = e+ cos 2</>+ e x sin 2tfi 
and e r x = — e + sin 2</> + e x cos 2</>. 



One can then construct three ellipticity CF's, C\(9) = ( e+(#o)e+(0o + 0) 



along the line connecting 2 points, t\ is invariant but e x changes sign, so C3 must 
always be zero. 

These three CF's have been measured for sources in the FIRST catalog that 
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are not 'point' sources, i.e for all sources that have a > 2". With this cut, there are 
about 60 sources deg -2 . The ellipticities used were those obtained before beam 
deconvolution, making the measured CF's only rough estimates of the true CF's. 
To obtain an estimate of the noise, the CF's were also measured for a randomized 
catalog in which sources are assigned random ellipticities. We found that the Ci's 
and C2's measured in the catalog were well above the noise. The GVs we measured 
were consistent with zero. While these results are encouraging, systematic effects 
are likely to contribute significantly to the signal and we are still studying these 
in detail. 

We have also investigated what CDM models predict for C\ and C^. A COBE 
normalised model with Q = 1 and h = 0.5 predicts values which are above the 
noise. Normalising to popular values of o~%, however, decreases the predicted signal. 
Further work is required before we will know if a weak lensing signal is detectable 
but if it were, it would provide a unique opportunity to probe mass fluctuations 
on scales that are very difficult to probe with optical surveys. 
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